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I. INTRODUCTION 



Compton scattering on the proton (7p) and the deuteron (7D) provides a clean and unique 
probe of nucleon electromagnetic structure, revealing information different to that obtained 
in electron scattering. During the past decade a number of excellent experimental programs 
have been dedicated to these two processes (see Refs. [1-5] and [6-10], respectively). At 
low photon energies, these experiments probe the static properties of the nucleon, such 
as its electric charge, magnetic moment, and polarizabilities. Above the pion-production 
threshold, the process becomes dominated by the excitation of resonances, most prominently 
the A(1232)-isobar. Many theoretical methods aim at understanding this process in both 
the low-energy and the resonance region. In particular, significant progress has been made 
recently using dispersion relations [11, 12] and effective Lagrangian models [13-15]. On the 
other hand, previous calculations using chiral perturbation theory (xPT) appear to work 
only at low photon energies — energies at or below the pion-production threshold [16, 17]. 
This present study attempts to extend these calculations above the pion threshold and 
into the Delta-resonance region. 

In the low-energy regime %PT seems to work extremely well. At next-to-leading order 
(NLO), i.e., third order in small momenta [= 0(g 3 )], heavy-baryon (HB) %PT predicts, for 
the electric and magnetic polarizabilities [18, 19]: 
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— -^teJ; 122xir4fm - (1) 

where a = e 2 /An ~ 1/137, $a — 1-26, f n ~ 93 MeV, m n ~ 139 MeV [47]. Since there are no 
Compton counter-terms present at 0(q 3 ), this is a genuine prediction of %PT. A prediction 
which — at least for the proton — is in remarkable agreement with recent extractions of these 
quantities from low-energy data, e.g. [20]: 

a p = (12.1 ± 1.1 ± 0.5) x 10~ 4 fm 3 , (3) 
P P = (3.2 ± 1.1 ± 0.1) x 10~ 4 fm 3 , (4) 

where the first error is statistical and the second represents the theory error of the fit to 
data. 

However, the agreement of the NLO HB^PT prediction with the experimental 7p cross- 
section data is good only up to photon energies 00 ~ 100 MeV [16]. The recent NNLO [0(g 4 )] 
calculation [17] agrees with experiment to slightly higher energies, but above uo ~ 120 MeV 
significant discrepancies begin to appear, most notably at backward angles. This is perhaps 
not surprising, since the Delta-isobar excitation is not included explicitly in this chiral 
expansion. And, as we shall argue, the breakdown scale of xPT without an explicit Delta 
is set essentially by the Delta-nucleon mass difference: 

A = Ma — Mtv ~ 293 MeV . (5) 

Thus, to extend the region of xPT applicability to uj ~ A, the Delta must be included as 
an explicit degree of freedom. 
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The Delta contribution for the Compton amplitude had already been analyzed using chiral 
effective Lagrangians with explicit Deltas in Refs. [21-23]. These studies focused mainly on 
nucleon polarizabilities. The predictions made in Refs. [19, 21, 22] are obscured by off-shell 
ambiguities, in particular by the so-called "off-shell parameters" which control the infamous 
spin-1/2 sector of the spin-3/2 Delta field. In a "reasonable" range for these parameters the 
Delta contribution to (3^ varies between and 14 x 10~ 4 fm 3 [22]. In contrast, Hemmert 
et. al. [23], to next-to-leading order in their small scale expansion (SSE) [24], find a result 
which is independent of the off-shell parameters, and thus is apparently a reliable prediction. 
But this prediction for the Delta-contribution to the magnetic polarizability is ^9x 
10 -4 fm 3 , in dramatic contradiction with experiment [25]. (For an attempt to remedy this 
using a "modified SSE" see Ref. [26].) 

In this work we include the Delta in the chiral Lagrangian in a fashion somewhat different 
to this literature. First of all, the Lagrangian is written such that the unphysical spin-1/2 
components of the Delta field decouple from observables [27, 28], hence no "off-shell param- 
eters" appear. This feature of our Lagrangian, besides removing the redundant parameters, 
allows us to dress the Delta-pole contribution in a manifestly covariant way. 

Furthermore, we set up our power-counting scheme so that it is both closely connected 
to the usual %PT without explicit Deltas in the low-energy region uj ~ m n , and extends to 
the Delta region uj ~ A. This is achieved by recognizing the hierarchy of scales: 

m w < A < A ~ 1 GeV , (6) 

where A stands for the "high-energy scale", the breakdown scale of the theory. Therefore, 
our scheme is rather different from the SSE of Refs. [23, 24] (see also Ref. [29]), where the 
Delta-nucleon mass-difference is assumed to be of order m n (i.e., A ~ m n ). A more detailed 
comparison of our scheme and the SSE will be given below. 

With the three-scale hierarchy (6) one in principle has two small expansion parameters: 
m^/A and A/A. We regard both of them as roughly the same size and so introduce a single 
small parameter: 

A m n 

5= A~X' (7) 
Note that this implies that m n scales as 5 2 . 

The validity of the scale hierarchy (6) and the expansion in powers of 8 (which we shall 
refer to as the 5-expansion) is to be judged by the success of the resultant effective-field- 
theory (EFT) description of processes involving the excitation of Delta. We regard the 
results we shall present here for 7p scattering as significant evidence in favor of this EFT 
expansion. 

To obtain the NLO result for 7N scattering in our scheme for both the low-energy and 
the Delta regions, the Delta-pole contribution to this process must be dressed, and then 
added to the NLO HB%PT result. This introduces two free parameters which characterize 
the strength of the 7 AT — > A transition: qm and Qe- Adjusting these parameters we find 
very good agreement with the experimental 7p differential cross section up to uj ~ 300 MeV, 
thereby extending the domain of applicability of chiral EFT into the Delta region. At the 
same time we also find reasonable values for the nucleon polarizabilities. 

In the next section we introduce the Lagrangian for the Delta and discuss its properties. 
Section III then describes the 5-expansion for Compton scattering on the proton. In par- 
ticular, we show that for uj ~ m w the power counting is very similar to that of HB;\PT, 
while for uj ~ A the power counting mandates resummation of the A propagator, thereby 
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dressing the A and giving it a finite width. In Section IV we summarize the elements of our 
calculation, and then in Section V we present and discuss the results of our NLO calculation 
for the differential cross section, as well as for the spin-independent polarizabilities a p , f3 p . 
We conclude in Section VI. 



II. THE CHIRAL LAGRANGIAN 

The pion-nucleon sector of the HB^PT Lagrangian is well discussed in the literature, see 
e.g. [19]. The terms relevant for our purposes are [48] : 



\Df7i a D» bc Tt c -\mln 2 + N 
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F^tt 3 + ..., (8) 



where 7r a represents the isovector pseudoscalar pion field, N is the isodoublet spinor field 
of the nucleon, r a are the isospin Pauli matrices, = <9 M — ieQA^, with Q representing 
the electric charge isospin operator (Q n = —ie ab3 , Qn = ^(1 + t 3 )), the electromagnetic 
field, F^ v = d ll A v — d v A^ = ^e^" 13 F a p, and k is the anomalous magnetic moment of 
the nucleon (k p ~ 1.79, n n ~ —1.91). 

Next we specify the terms involving the Delta field. Describing the Delta field by an 
isospin-3/2 spin-3/2 Rarita-Schwinger (RS) vector-spinor A M (x), we write the Delta piece 
of the chiral Lagrangian in the following form: 

C = Crs + £?riVA + ^7 AT A + • • • , (9) 

C RS = A, d a - M A r v ) K , (10) 

£.na = 7^rNTlr vX {dAu) dxn a + H.c. , (11) 

*J-k M A 

C ^ = 2M N[ M N + M~ - ^ F> ") ^ + H ' C ' < 12 » 

These are the free spin-3/2 Lagrangian, and the nNA, ^NA couplings respectively. Here 
T a are the isospin 1/2 to 3/2 transition matrices satisfying TjTf, = ^5 a b — \ie a b c T c . 

We have kept only the couplings that are linear in the A field and lowest order in the 
pion and the photon fields. In principle, there are many other couplings (ttttNN, ttttNA, 
7AA, etc.), represented in Eqs. (8) and (9) by the dots, that are required by the chiral and 
electromagnetic-gauge symmetries. However, they are not relevant for our calculation at the 
order we will be considering. 

For the purposes of power counting we rearrange the interaction Lagrangian according 
to the number of small quantities (momentum, pion mass, and factors of e) that each term 
carries: 

Ci = ££ w , 

i 

£W = -^N^ ■ dit a )T al5 N + *^^N 1 ■ Aitt 15 N + eNQ Nl ■ AN + C wNA , 
£(2) = ^N^NF^ + ±(ieQ n nA-dn + E.c.) + e 2 QlA 2 nn + d 9 ™l, (13) 
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r(3) — r^ 9E "> 

P 2 

£( 4 ) = p F^ir* 



A. Spin-3/2 gauge invariance 

It is important to note that our iVA couplings, besides being chiral and gauge invariant, 
are invariant under the following local (gauge) transformation of the spin-3/2 field: 

A^x)^A^x)+d^e(x), (14) 

where e is a spinor. This invariance ensures that the spin-3/2 field has the correct number 
of spin degrees of freedom (i.e., 2s + 1 = 4), cf. [27, 28]. 

As the result of this additional symmetry, any vertex involving a Delta field, r^(p, . . .), 
with n being the vector index and p the 4-momentum of the Delta, will obey the transversality 
condition: 

p M r"(p,...) = o. (15) 

Using the well-known form of the spin-3/2 propagator: 



S^(p) = 1 



f 1 — M{y 



(16) 



it is easy to show that the spin- 1/2 sector of the RS propagator decouples [28] and one may 
equivalently use the following propagator: 

~ SM = J^MfV 2,(p) (17) 

where VfJ 2) {p) = - §7^ - 3^2(^7^ + VyTfvf) is the spin-3/2 projection operator. 
As a matter of fact, it is then also possible to replace the vertices as follows: 

r l >(p,...) = v$, m r»(p,...). (is) 

In this theory P and r are completely equivalent. Nevertheless, vertices P are sometimes 
more convenient in actual calculations. For example, the irNA vertex from Eq. (11): 

r»(p,k) = (g/M A )^Pakp, (19) 

where g = tiA/Zfn and k is the pion 4-momentum, can be replaced by 

r,(p,k) = (g/M A )p'V% 2 (p)k», (20) 

where we have used: p^V% 2 (p) = = ^Vf/J 1 (p) . Furthermore, using 

V 3/2 V 3/2 = V 3/2 and 

\]/>,V 3 / 2 (p)] = 0, Delta-exchange amplitudes are computed effortlessly, e.g., 

r"(p, k') s»„{p) r»( P , k) = r"(p, k>) ~s, v { P ) r"(p, k) = j^^ffJ 2 \p) W ■ (21) 
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B. Relation to conventional A couplings 



Our iVA couplings are rather different from the usual ones of, e.g. Refs. [13, 19, 22-24]. 
As a rule, standard couplings do not have the spin-3/2 gauge symmetry (14). Exceptions 
are the 7AA coupling of Jones and Scadron [30], which obviously satisfies (15), and the 
couplings used by Kondratyuk and Scholten [15]. We have adopted the Jones and Scadron 
convention for the magnetic (Ml) coupling, gM in Eq. (12). 

Other conventional couplings, including the Gi, G2 representation of the 7AA vertex, do 
not have the spin-3/2 gauge symmetry. As a result they involve the unphysical lower-spin 
sectors of the spin-3/2 field, and hence observables become dependent on arbitrary "spin- 
1/2 backgrounds" associated with "off-shell parameters" of the Delta. Other pathologies 
(see Ref. [27] and references therein) — all of which can be traced back to the fact that the 
couplings violate the degrees-of-freedom-counting constraints of the free theory — also occur 
in these theories. 

One can establish a relation between the "inconsistent" and "consistent" couplings using 
field transformations [31], but this relation holds only in perturbation theory [49] and so is 
not strictly applicable when resummation of the Delta contributions is necessary, as is the 
case in the computation of the next section. 

III. COMPTON AMPLITUDE IN THE 5 EXPANSION 

In Compton scattering the momenta of the particles are characterized by the photon 
energy uo. For very low photon energies pions can be "integrated out" of the theory with 
all non-analytic effects associated with their production being replaced by a power series in 
uo/m^ (see, for instance Refs. [32]). Clearly, the condition for this EFT to be effective is 
uo <C m n . If, instead, we want to develop an EFT for uo ~ m n we must treat both uo and m n 
as low-energy scales, which means that pions must appear in the theory as explicit degrees 
of freedom. 

Similarly for the next relevant scale: A = Ma — 293 MeV. A theory which does 

not treat it as a low-energy scale is effective only for uo <C A. %PT without explicit Deltas 
is an example of such an EFT. There only is treated as a small scale, and it is assumed 
that uo ~ m v <C A. To extend this description to the Delta region, uo ~ A, we need to adopt 
A as a low-energy scale and include the Delta as an explicit degree of freedom. Thus, we 
naturally arrive at the scale hierarchy: 



This hierarchy complies with the assumption of %PT: m n <C A, and so still gives the 
dominant effects in the theory if uo ~ m^. At the same time (22) allows us to extend %PT 
to the Delta region. 

In developing our power counting below, we will often keep the dependence of amplitudes 
on m n and A explicit, so that the behavior of the amplitudes in the (independent) chiral 
(m^ — > 0) or the large- N c (A — > 0) limits is manifest [33]. Nevertheless, for the purposes of 
assigning an overall size to the amplitude arising from a particular graph or set of graphs 
we would like to have one expansion parameter: 



m n < A < A . 



(22) 




A 



(23) 
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where we conservatively adopt A pa 600 MeV, the scale introduced by the excitation energy 
of the next baryon resonance. In fact, A will represent not only this scale but all of the 
various high-energy scales, such as m p , Mat, Ma, and 47r/ 7r . Obviously, in this counting, A 
scales as 5, while scales as 5 2 . 

While 5 is of order one half, the expansion in powers of 5, is, in principle, no worse than 
%PT, which is an expansion in m n /A, or the SSE [24], which is an expansion in powers of 
A/A. Note that, since Eq. (23) is not necessarily true in worlds with other values of N c , m q , 
etc., once that equation is employed, the connection to the limits m q — > and N c — > oo is 
lost unless the chiral and large- N c limits are taken simultaneously with m^N 2 held fixed. 

We assign to each graph an overall 5-counting index, a, which simply tells us that the 
graph is of size e 2 5 a /A. Because we deal with two different low-energy scales in our EFT, 
the index a has two different expressions, depending on whether the photon energy uj is 
in the vicinity of m n or A. For a graph with L loops, N n pion propagators, N N nucleon 
propagators, N A Delta propagators, and Vi vertices of dimension i, the index is 

a= (2a xPT -N A , uj~m; 
{ a x PT ~ n a , uj ~ A, 

where q; x pt = Hi iV% — 2 + 4L — N N — 2N n is the index of the graph in xPT without explicit 
Deltas. 

In deriving this power-counting we have used the fact that no graphs containing vertices 
with powers of m 2 or A occur up to the order to which we work. Such vertices do arise in 
higher-order graphs though. In general then, a vertex with j derivatives, k insertions of the 
quark mass, and I powers of the A-N mass difference scales as u; J m 2 ^ A', and so has overall 
dimension i = j + 2k + I. Denoting the number of such vertices by VjM, the 5- index of an 
arbitrary graph is: 



a = < 



f E(2j + 4A; + l)V jkl + 2{AL - N N - 2N W - 2) - N A , uo ~ m w ; 

jkl 

E(j + 4A; + l)V jkl + AL-2-N N -2N w -N A , oo ~ A, 



which obviously reduces to Eq. (24) if only vertices with k — I — are present. 

In the region uj ~ A, there is an important exception to this scaling rule: graphs that 
are one-Delta-reducible (OAR), such as those in Fig. 1, scale not as 5 a but as 

NoAR 

6 a ■ 



/ 1 \ "OAR 



where N A in the equation for a now counts only the one-Delta-irreducible propagators, while 
Noar is the number of OAR propagators. In the low-energy region this does not affect the 
power-counting, however in the region uj ~ A these graphs can be dramatically enhanced. 
This forces us to resum all the OAR contributions, which amounts to dressing the Delta 
propagator, thus ameliorating the divergence which otherwise occurs at the Delta pole, and 
producing a width for the Delta of roughly the experimentally-observed size. 

Details of the dressing are given below. By definition, once dressing is performed a OAR 
graph can have only one Delta propagator, and such a graph then scales as: 



*" • < 26) 
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(a) (b) (c) (d) 

ki Lj U 

(e) (/) (s) 

FIG. 1: Eight graphs which are one-Delta-reducible (OAR) and so become enhanced in the region 
lo ~ A. The sliced vertex is the Ml jNA vertex, while the sliced and diced 7./V A vertex is the E2 
coupling from C^ 3 \ Solid dots represent couplings from 



where E is the self-energy. The expansion for E begins at 5 3 , and so in the domain \u — A| ~ 
5 3 the OAR graphs are enhanced by 5~ 2 over the value expected from Eq. (24). Thus, the 
correct index of a OAR graph in the region uj ~ A is 



a = a 



2. 



(27) 



As a result, for instance, the s-channel-pole Delta graph of Fig. 1(a), which is the simplest 
OAR graph, is promoted from NNLO in the low-energy region to LO in the Delta region. 

The rest of this section is organized as follows: before giving a detailed explanation of 
Compton counting in the (^-expansion we make a few comments on how our scheme compares 
to standard HB^PT and to the SSE of Ref. [23, 24]. We then discuss power counting for 
the low-energy region uj ~ m n . In subsection III C we explain the central issue for the 
higher-energy domain uj ~ A: the dressing of the Delta pole. Then in subsection III D we 
elucidate the impact of this dressing on the counting for Compton scattering graphs. 



A. Comparison with HB^PT/SSE 

In HB^PT the Delta is not included as an explicit "low-energy" degree of freedom in the 
Lagrangian. Instead it is integrated out of the theory, producing a low-energy theory that, 
in principle, breaks down for uj ~ A. Power counting of graphs is then performed in terms 
of the index q, where 

* = X~X' (28) 

where A is usually assumed to be of order 1 GeV, although the omission of explicit Deltas 
suggests instead A ~ A. 

Hemmert et al. [23, 24] introduced the "small-scale expansion" (SSE) (see also the earlier 
Ref. [29]), where the EFT expansion parameter is: 

The SSE treats and A as the same scale, and hence the Delta must be included explicitly 
in both energy domains: uj ~ m w and uj ~ A. 
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This overemphasizes the importance of the Delta somewhat at low energies. In contrast, 
in the low-energy region, the ^-expansion amplitude is akin to that of HB^PT. In the region 
uj ~ A the dressing of the Delta implemented here is not performed in either HB^PT — 
naturally, since Deltas are "high-energy" degrees of freedom — or the SSE — as for uj ~ A ~ 
all ttN loop effects are a small correction to the "bare" Delta propagator. 

The table below summarizes the relationship of the ^-expansion to HB%PT and the SSE. 



Expansion 




A/ A 


HBxPT 


q 


1 


SSE 


e 


e 


(5-expansion 


5 2 


5 



TABLE I: The three different expansion discussed in the text. In all three cases the small expansion 
parameter is of order 1/2, and A is the breakdown scale of the theory. 



B. Power counting for u ~ 

Here we make the identification uj,m n ~ 5 2 . Graphs without Delta propagators then 
scale exactly as in x?T, but with small momenta q = 5 2 [50]. The general index of such 
a graph is then given by Eq. (24). The leading contribution then comes from the sum of 
the relativistic nucleon Born graph (with V\ = 2, Njy — 1, L — N n = = 0) depicted in 
Fig. 2(a) and its crossed partner. Both graphs behave as e 2 jui ~ 5~ 2 as uj — > 0. But the 
divergent parts cancel in the sum, as the low-energy theorem tells us they must [34]. The 
dominant term for small uj is given by the "Thomson amplitude" : 

T ™ = < 3 °> 

with e and e' the photon's initial- and final-state polarization vectors. This, obviously, is 
0(5°). 



y y y 




(a) (b) (c) (d) 

FIG. 2: The four relativistic tree-level graphs without Deltas that are included in our calculation. 
(We also include graphs generated from these graphs by crossing and/or time reversal.) The dot is 
the leading-order jNN coupling, while the sliced vertex is the anomalous-magnetic-moment vertex 
from C (2 \ The square indicates the ir° —> 77 vertex from which generates the chiral anomaly. 

When the expansion of the relativistic graphs Fig. 2(a) and 2(b) in powers of uj is made 
there are also pieces ~ e 2 uo = 0(5 2 ). These form part of the NLO amplitude. The rest of 
the NLO amplitude is obtained from graphs which have index a = 2: nucleon tree graphs 
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with the anomalous magnetic moment coupling (i.e., V 2 = 2, Nn — 1, L — N n = = 0), 
see Fig. 2(c), and the ir° exchange graph 2(d), involving the WZW anomaly, which has 
Vi = V 4 = 1, N n = 1, L = N N = N A = 0. 

Next we consider 7rN-loop contributions to 7N scattering. After making a heavy-baryon 
expansion of our relativistic Lagrangian, in order to avoid difficulties with the appearance 
of the scale inside loops [51], we construct the leading loop graphs from vertices in 
This yields the graphs of Fig. 3, together with their crossed partners, as reviewed in 
Ref. [19]. These graphs are specified by L — 1, N N = 1, = and either V\ — V 2 — 2, 
N v = 3 or Vi = 2, N v = 1 or V x = 2, V 2 = 1, N w = 2, and hence all have a = 2. They are 
the only loop graphs with this counting index if we adopt Coulomb gauge and employ the 
heavy-baryon expansion. Explicit computation [17, 19] reveals that the sum of these graphs 
indeed produces a Compton amplitude that behaves as: 

2 2 

T (*N loop) = e U F (D ( M ? ( 31 ) 

where is a non-analytic function whose form is given in detail for the various possible 
spin and polarization structures in Refs. [17, 19] and in Appendix A. has the property 
that F^ ~ 1 for uo < m f . A crucial feature of Eq. (31) is the fact that the sum of these 
leading loop graphs is proportional to uo 2 . This is also a consequence of the low-energy 
theorem [34]. 




(») (<=) M 

FIG. 3: The dominant 7rN-loop contributions to 7N scattering (crossed and time-reversed partners 
are not shown, but are included). 

The counting formula (24) indicates that loop graphs with insertions from the second- 
order xPT Lagrangian 

£(2) 

are down by two further powers of 5, being of 0(<5 4 ). Relativistic 
corrections to Eq. (31) are suppressed by u/M, and so are also 0(5 4 ). Some loop graphs 
at 0(5 A ) require renormalization, and the corresponding counterterms must be included. 
Meanwhile, graphs with two 7rN loops are 0(<5 6 ) in this counting. Thus — at least in this 
energy domain — it is not until 0(5 6 ) that two-pion intermediate states contribute to the 7N 
amplitude. And graphs involving additional 7rN rescatterings are similarly suppressed. Con- 
sidering more loops and/or insertions with more derivatives only serves to further increase 
the 5-index of graphs. Thus unitarity (in both the 7rN and 7T7tN channels) is violated in our 
calculation, but the violation is always an effect of an order beyond that at which we work. 

Graphs containing the Delta begin to contribute at 0(5 3 ). The tree graph with two Ml 
"fNA vertices — see Fig. 1(a) — has a = 3 (V2 = 2, ]Va — 1, L — Njy = N n ). Meanwhile, the 
counting of the one-7rA-loop graphs, Fig. 4, is analogous to that for the nN loops, the only 
difference being that now Nn = 0, TVa = 1 instead of = 1, A^a = 0. This results in the 
sum of the graphs in Fig. 4 being 0(<5 3 ), i.e. scaling as: 

T (An loop) _ e2 (1) ( 
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where is a non-analytic function which is of order 1 for m,/ A ~ 5, and u/A<l. Equa- 
tion (32) is consistent with the low-energy theorem, and agrees with the explicit computation 
performed for these loops in Ref. [24]. 




w w 



(6) (c) (d) 

FIG. 4: The dominant 7rA-loop contributions to 7N scattering. Again, graphs generated from 
these by crossing and/or time-reversal are not shown. 

In summary: in the ^-expansion the Thomson term is the leading mechanism for Compton 
scattering on the nucleon at low energies, u ~ m n . In this region pion loops are suppressed 
by one power of u ~ m n ~ 5 2 , exactly as in HB^PT. If explicit Deltas are included in the 
theory the leading Delta-pole and A7r-loop graphs are suppressed by 5 3 relative to leading. 
They thus occur one order higher in the (^-expansion than the Nn loop graphs of Fig. 3. 
They are, however, still one power of 5' 1 larger than graphs arising from insertions in 
%N loop graphs. 



C. Dressing the Delta 

The key issue for the theory in the region u ~ A is the treatment of the Delta pole. 
One-A-reducible (OAR) graphs must be resummed in order to remove the divergence which 
otherwise occurs when j> = 

Formally, all OAR graphs can be summed via the series: 

SM = S$(p) + S%(p) E^'{p) S$l(p) + ..., (33) 

where S^ u (p) is the full one-Delta- irreducible (OAI) Delta self-energy, and S fiu (p) (Sj®J(p)) 
is the dressed (bare) Delta propagator. 

The function has a (^-expansion of its own: 

^ = + Eg + • • • • (34) 

This expansion begins at 0(5 3 ), with the graphs depicted in Fig. 5, together with the 
counterterms necessary for their renormalization. Insertions from C^ N generate effects in 
£( 4 \ These effects include relativistic corrections to the leading heavy-baryon result S^ 3 ^. 
Two-loop contributions to the self-energy — including the leading effect of the Nun channel — 
first occur in £( 5 ), and are thus smaller by 5 2 than the dominant piece of S. 
If: 

|cu-A|~<5 3 (35) 

and we keep only the third-order piece of the self-energy then all terms on the right-hand 
side of Eq. (33) are the same order. A perturbative expansion of the right-hand side — which 
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(a) (b) 

FIG. 5: 7rN and ttA contributions to the A self-energy S. The vertices are from and so both 
graphs are 0(5 3 ). 

is certainly valid for \uj — A| ~ 5 2 or larger — is no longer appropriate. Instead, if (35) holds, 
the whole series must be resummed, giving: 

Eq. (35) then defines precisely what we mean by u ~ A. 

In fact the most general Lorentz-covariant form of S tJ,u (p) is rather complicated: it 
contains up to 10 independent scalar functions. As a result the dressed propagator does 
not generally have the form (36). This is a consequence of using "inconsistent" spin-3/2 
couplings — ones that do not obey the symmetry under (14). If, however, couplings which 
are consistent in that sense are used then the Delta self-energy can be written as: 

27"" (p) = S{p)V {3/2) ^{p), (37) 

with 27(p) akin to the usual fermion self-energy: 27(p) = A(p 2 )^> + B(p 2 ), where A and B 
are scalar functions. Dressing then affects only the spin-3/2 piece of the propagator and the 
form (36) results. The divergence at p 1 = M A is ameliorated, and no further resummation 
is necessary. 5-counting indicates that the effects of £( n ) for n > 4 can be included by 
perturbing around the propagator (36). 

Delta propagators of this form have been used in other authors' extensions of chiral 
perturbation theory to the resonance region [36-38], although in these works it is not clear 
why only the spin-3/2 sector is dressed. Note that in contrast to the work of, for instance, 
Ref. [37], we do not dress the nucleon pole by nN loops. Arguments analogous to those of 
this subsection suggest that nucleon dressing is only necessary from a power-counting point 
of view for u ~ 0, and there S(p) is purely real. As we shall now see, after renormalization 
the real part of 27^ 3 ^ plays a negligible role in the propagator (36). 

In Eq. (36) the quantity indicates that we are resumming the renormalized third-order 
Delta self-energy. The explicit renormalization of this quantity will be performed elsewhere. 
Here we make a more general argument which constrains the form and importance of any 
renormalized self-energy appearing in Eq. (36). 

First, observe that the general Lorentz structure of the self-energy £ results in: 

s ^ f) - -j^m?-/^ (38 > 

with Z and M scalar functions of p 2 . After mass, wave function, and coupling constant 
renormalization these can be written as: 

Z(p 2 ) = 1 + (p 2 - M 2 A )f z (p 2 )+iImZ(p 2 ); (39) 
M(p 2 ) = M A + (p 2 - M 2 A ) 2 f M (p 2 ) + i Im M(p 2 ), (40) 
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with fz and /m real functions of p 2 . 

Substituting these forms into Eq. (38) we find that: 

S ( V ) = l + ^m^(p 2 ) p W2 )() + (l\ m 

^ v[P) j> — Ma — % Im M(p 2 ) ^ [P) + U \a) ' 1 j 

In the region (35) the pole piece is 0(5~ 3 ), while the remaining "background" terms are 
0(5°). Thus any corrections to S^ip 2 ) from ReZ(p 2 ) or ReM(p 2 ) are three powers of 5 
beyond leading. Corrections from lmZ(p 2 ) are equally suppressed, since it too is 0(5 3 ). 

Thus, up to corrections which are N 3 LO, it is sufficient to compute only Im M{p 2 ). Pro- 
vided that uo < A + mn this comes exclusively from Fig. 5(a). A straightforward calculation 
gives: 

where k is the on-shell value of the pion three-momentum: 

k = {[s- (M N + m n ) 2 ][s - (M N - m w ) 2 ]/(4s)} 1/2 ~ 5. (43) 

Thus, the width is 0(5 3 ), as promised. 

The final form of the resummed A propagator is then: 

If this propagator appears in a OAR 7N graph and (35) is satisfied then it scales as 5~ 3 . 



D. Power counting for uj ~ A 

The effect of this modified scaling for the Delta propagator is that (dressed) OAR graphs 
become the dominant effects for uo ~ A. Their 5-index is given by Eq. (27). The Delta-pole 
graph Fig. 1(a), with Ml 7AA vertices (i.e., Vi = 2, AA = 1), has a = —1 in the region 
(35) and gives the leading contribution there. The graph of Fig. 1(h) with one E2-coupling 
(i.e., V2 = 1, V3 = 1) has a — 0, and hence is of next-to-leading order (NLO) if uj ~ A. 

Meanwhile the OAR graphs of Fig. l(b)-(g) (and their time-reversed partners) are char- 
acterized by L = N N = N A = 1 and either V x = V 2 = N n = 2 or V 2 = N n = 1, V x = 2. 
They too have a = and contribute at NLO. These loop graphs are divergent and must be 
renormalized. This is achieved via graph 1(h). The loop effects may then be included in the 
calculation by the use of an energy-dependent E2-coupling: g E — > Qe^s). The leading effect 
here again arises from the imaginary part of the loops, hence: 

( , , . ( 9Ah A \ s + M 2 N -ml M N k 2 2M N {M N + M A ) ^ 

(45) 



where w 7 = (s — M^ r )/2 v /s, uo k = \jm\ + k 2 , while k is given by Eq. (43), and Qi is the Ith 
Legendre function of the second kind. 

NLO effects can also be obtained by considering corrections from to the leading Delta 
self-energy X®. More complicated electromagnetic couplings and higher-order terms in E 
lead to effects of NNLO in the (^-expansion of OAR graphs. 
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What is the 5-index of graphs without a Delta pole? OAI graphs such as Aty and Ntt 
loops obey Eq. (24) — they are not enhanced. They retain a positive S index, and so are at 
least NNLO in this counting. 

One might wonder how to reconcile this with Eq. (31) which seems to suggest that the 
^-scaling of the dominant tcN loops will be, for uo ~ A: 



, > 2 A 2 

2 w 2 L -^ 2 rO 

e — = e — = e 5". 

VTLjr A A 



(46) 



This conclusion is erroneous because is not 0(1) if uo/m^ is large. If uo/m^ ^> 1 the 
loop functions should be expanded about the large-u; limit, not the small uo one, and doing 
so results in: 



rp{jrN loop) 



(4tt 2 /0 



2 4 
777, YYl 

CiUJ + dim^ + C 3 — - + + 



UO 



UO" 



(47) 



Assigning uo ~ A, m 7r ~ A 2 we see that this series, rather than one in increasing powers of uo 
is the correct one for the "medium-energy" regime uo ~ A ^> m n . This is completely opposite 
to a polarizability expansion in increasing powers of uo. In our approach uo ~ A is sufficiently 
far from threshold that a power-series expansion around uo = of Delta contributions is not 
very enlightening. 



IV. IN PRACTICE 

A. Defining the NLO calculation 

To perform a complete NLO calculation in the whole energy region, < uo < A, we 
include all of the nucleon pole and nN loop graphs of Figs. 2 and 3, together with their 
crossed partners. To these we add the OAR graphs of Fig. 1, which in the low-energy region 
contribute at NNLO and above, but for \uo — A\ ~ 5 3 give effects of leading and next-to- 
leading order. We also include the OAR graph with two g E 7AA couplings, even though it 
is formally NNLO. 

We keep all of these effects in both kinematic regions. Note that this means we are 
always keeping contributions that are, strictly speaking, beyond the order to which we 
work. This is done in order to provide a smooth transition between the two different photon 
energy domains. For the same reason, in both regions we always use the resummed Delta 
propagator (44). 

The power counting of Sec. Ill C indicates that at NLO we must include effects due to the 
0(S 4 ) piece of S, S^. The heavy-baryon graphs which contribute to the 0(e 4 ) self-energy 
in a heavy-baryon calculation with explicit Deltas [40] are shown in Fig. 6. The relativistic 
calculation of the Delta width that led to Eq. (42) already includes the effects of the graphs 
(a) and (b). As for graph (c), this graph gives a contribution to £ which behaves as [40]: 

The imaginary part of this graph is proportional to uokm^ (with k given by Eq. (43)), so 
while it is 0(e 4 ) in the SSE, in our counting it is 0(<5 6 ), and so well beyond the order to 
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(a) (b) (c) 

FIG. 6: Heavy-baryon expansion graphs that represent 0(e 4 ) 7rN-loop contributions to the A 
self-energy S. The sliced squares are vertices from the second-order ttNA HB Lagrangian whose 
coefficient is fixed by Lorentz invariance. The sliced circle represents a vertex from C^na with a 
coefficient which is a priori undetermined (see Ref. [40]). 

which we work. Thus the result (42) is is already accurate up to corrections of relative order 
5 2 . 

The w-channel Delta-pole graph is NNLO throughout this entire energy region and is not 
included in our Compton amplitude. Therefore, by looking at its effect on cross sections we 
can estimate the importance of NNLO contributions. 



B. Covariant decomposition of the Compton amplitude 

To compute the amplitude of Compton-scattering on a spin-1/2 target in a manifestly 
Lorentz and gauge invariant form, we specify it in terms of eight scalar amplitudes A-y . . . A 8 
(instead of the usual six): 

M fi = e 2 u(p') j2M*,t)OrtfW)£v(q)u{p) (49) 
i=i 

where p', p (and q', q) are the final and initial 4-momenta of the nucleon (and photon), 
respectively; u is the free nucleon spinor, is a modified photon polarization vector: 

P ■ e 

with P = p + p' . Tensors Oi are given by: 



or 


= -9^ 


Of 


= q"q' U 


or 


= ~y u 


Of 




or 


= q"q' a i au - -rw 


or 


= qW ~ l" a q a q' 


or 


= q"q V (q'-l- q) 


or 


= il^tM ■ 



(51) 



Mandelstam variables s, t, u are denned as usual: 

s = (p + qf = M 2 N + 2p-q 

t = (q-q) 2 = -2q-q (52) 
u = ( p - q 'y = M 2 N -2p-q' 
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where we have used the onshell conditions: q 2 = q' 2 = 0, p 2 = p' 2 = Mjy. Note that 
P.q = l(s-u) = P-q'. 

The representation (49) is obtained by writing down the most general covariant structure 
(for the onshell situation) and imposing the electromagnetic current-conservation condition. 
Thus this representation incorporates both covariance and gauge-invariance in a manifest 
way. 

The amplitudes A* are most easily computed in the following Lorenz-invariant gauge: 

P . £ = = P ■ e' . (53) 
This condition can also be achieved in the Coloumb gauge (e = = e' ) by going to the 

— * 

Breit frame: P = 0. 

In the Coloumb gauge, the structures 0±-Oq exactly match the ones in the standard 
decomposition (e.g. [19, 34]), while O-j and Og can be reduced to linear combinations of O3, 
O5 and Oq. 

7 = (uj 2 + it) [-t 3 + 5 ] - 0J 2 O 6 , (54) 

8 = ±u(-t0 3 + 5 - 6 ). (55) 

In particular, the results for the nucleon-exchange (Born) graphs, the anomaly graph, and 
the Delta exchange graphs are specified in terms of the amplitudes Ai, . . . , A 8 in Appendix 
A. As our calculation of these graphs is fully relativistic, it differs from that of HB^PT in 
Refs. [19, 23] by terms of 0(5). 

Meanwhile, for the 71N loop graphs depicted in Fig. 3, we have used the 0(q 3 ) HB^PT 
result of Refs. [17, 19]. This is done in order to avoid difficulties with the treatment of 
the nucleon-mass-scale Mn inside the loops. The amplitude obtained in a fully- relativistic 
calculation of these loops would differ from the result used here by terms of 0(u/M) and 
0(m n /M), i.e. terms down by 5 2 . The loop contributions are given in Appendix A. 

We do include one particular relativistic effect because we write the photon energy uj that 
appears in the loop functions as u = ^Js — Mat. The standard choices for uj — cm. photon 
energy [19], Breit-frame photon energy [17] — differ from this by terms which are of N 3 LO in 
our counting. We adopt this prescription in order to ensure that the Compton amplitude's 
7iN cut occurs at the correct value of s. 

V. RESULTS AND DISCUSSION 

A. Differential cross sections for 7p scattering 

There are no 7N contact terms at NLO in the (^-expansion. This leaves us with three 
EFT parameters which must be fixed using the data: qe, Qm, and Ha- For Ha we adopt the 
phenomenological value h A = 2.81 {fl N ^/^ — 0.35), which corresponds to a Delta width 
r(M\) ~ 111 MeV — consistent with the range given by the Particle Data Group [39]. This 
value of h-A is roughly 5% larger than that obtained from the large- N c relation Ha = ~j^9a- 

This leaves us with two free parameters: qm and qe of the 7A^A coupling (12). They 
represent the strength of, respectively, the Ml and E2 7A — > A transitions. 
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In principle these are relatively well known from pion photoproduction. In particular, 
their ratio is related to the Rem = E2/M1 ratio: 

Rem = 7: 777 \, v r- x 100 %. (56) 

2g M (M N + M A ) - g E A v ; 

The determination of this ratio has recently been the subject of experimental programs 
at JLab and MAMI. The present PDG value is R EM = (-2.5 ± 1.0)% [39]. However, 
this number is measured only indirectly through the extraction of the ratio of the pion- 
photoproduction multipoles at the Delta-resonance position. These multipoles are affected 
by a number of background processes and the relationship (56) is only strictly true at leading 
order in the ^-expansion for the 7P^ 7rp amplitude. To fully understand the constraint that 
pion photoproduction places on Qe/qm in this EFT a higher-order calculation of 7P^ 7rp 
using the ^-expansion is necessary. (For an attempt to compute this process in the SSE see 
Ref. [41].) Absent such a calculation, here we regard gu and g E as free parameters and fit 
them to get the best agreement with 7p cross-section data. An important future test of the 
usefulness of the (^-expansion will be whether the resultant value for g E jgu is ultimately 
consistent with that found from pion photoproduction data using the same framework. 

The results for the differential cross section at several different energies are presented in 
Figs. 7 and 8. The long-dashed orange curve represents a calculation which includes only the 
Born graphs of Fig. 2. The dashed blue curve is the result when the %N loops of Fig. 3 are 
added, and so gives the 0(q 3 ) prediction of HB%PT. Finally, the complete NLO calculation 
in the (^-expansion is represented by the solid red curve. The sharp rise at backward angles 
as uj increases past the pion-production threshold is now reproduced in the theory. This 
sharp rise is very difficult to obtain in x?T without explicit Deltas. 

In Figs. 7 and 8 we also show a theoretical error band, demarcated by the red dots. This 
is obtained by estimating the size of the NNLO contribution in the following way: 

T(theor. err.) = e> ■ e x { W ^ A ' w ~ ™* (57) 

v ' M N [uj, ~ A . 17 

The error band in Figs. 7 and 8 is plotted with the border between the two kinematic domains 
at uj — 200 MeV. The band is an estimate of how far we expect the NNLO corrections to 
change the results. It may overestimate the theoretical error since our calculation already 
includes a number of NNLO contributions: relativistic effects coming from tree graphs, 
Delta-exchange — which is NNLO for uj ~ — and irN loops — which are NNLO for uj ~ A. 
The fit of the cross section favors the following values for the 7^ — > A parameters: 

g M = 2.6 ±0.2, g E = -6.0 ±0.9, (58) 

The solid line in Figs. 7 and 8 gives the result for the central values of gu and g E . The 
uncertainty in Eq. (58) is found by varying these couplings until the experimental cross 
section in the Delta region is no longer within the theoretical error band. 

The resulting value of gu is consistent with the large N c value: gu = 3^(1 + K p) — 2.63, 
while the value of g E is considerably different from phenomeno logical expectations. For 
instance, from R E m = —2.5% and Eq. (56) one would expect g E = —1. This problem 
dramatizes the importance of performing an analogous calculation of pion photoproduction 
to check the consistency of the (^-expansion. 
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Next we would like to note that, although formally the imaginary and real parts of Qe{s) 
of Eq. (45) are the same order, numerically the imaginary part is smaller by at least a factor 
of six. ImgE(s) has a negligible impact on the angular distributions shown in Figs. 7 and 8. 

The NLO prediction for the differential cross section's energy-dependence is shown in 
Fig. 9 for a scattering angle of 90°. The solid purple line is the result when only A-pole 
mechanisms are included and qm and g% are chosen according to (58). The solid red line 
gives the result of our NLO calculation. Individual contributions from Born graphs and nN 
loops are given by the orange dashed and green short-dashed lines respectively. The 0(q 3 ) 
prediction of HB%PT is represented by the blue dashed line. This figure shows our NLO 
calculation of the Delta width is in good agreement with the data. This lends support to our 
adoption of the value Ha = 2.81. However, it must be pointed out that since our EFT was 
only designed for < uo < A the agreement at the higher energies which is seen in Fig. 9 is 
probably somewhat fortuitous. 

B. Polarizabilities 

The results for the nucleon polarizabilities are worked out in Appendix B. The leading 
contribution to polarizabilities is from ixN loops and hence is exactly the same as in HBxPT, 
see Eq. (1). On the other hand, our results for the Delta contributions differ from earlier 
ones in several ways: 

1. First of all, in contrast to the results of Refs. [13, 22], they are independent of "off-shell 
parameters" . 

2. The leading SSE result of Ref. [23] for the magnetic polarizability is also free of off- 
shell parameters (in the SSE they enter at NLO). Formally, the result of Appendix B 
agrees with that found in Ref. [23] — apart from a higher-order relativistic effect which 
is numerically only of order 10%: 

R (A) this work _ 2 ( A ) SSE 

Pn ~ 2 + A Pn 

However, our fit to the cross section prefers a significantly smaller value of the 7./VA 
coupling than was used in Ref. [23] and this leads to a markedly smaller numerical 
result for (3^\ 

3. Perhaps most importantly, our NLO calculation includes neither w-channel Delta- 
exchange nor 7rA loops, since both are NNLO in the 5-expansion. In fact, the Delta 
contributions quoted in Appendix B stem from both the s- and w-channel Delta ex- 
change graphs. However, the NLO calculation presented above has no w-channel graph. 
Thus to find the polarizabilities that correspond to the cross-section calculation of the 
previous subsection one must halve the Delta pieces of polarizabilities given in Ap- 
pendix B. 

All of these effects produce a Delta-contribution to the magnetic polarizability that is 
significantly smaller than that found in previous EFT calculations with explicit Deltas [21- 
23, 25]: 

^'- (M^yA -P- 7 ^*"^' (59) 
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The Delta contribution to the electric polarizability comes solely from the s-channel Delta- 
pole with two Qe couplings: 



af) = a9E = (-2.0 ±0.7) x 10- 4 fm 3 . (60) 

v (Mjv + Ma) 3 v 1 v 1 

The observant reader will have noticed that the inclusion of A contributions in the polar- 
izabilities is not, strictly speaking, consistent at NLO in our power counting. The s-channel 
Delta pole is, like its w-channel counterpart, NNLO for uj ~ 0. In spite of this we have in- 
cluded the Delta contributions in the results shown in Table II. The numbers there represent 
the polarizabilities corresponding to the cross-section calculation already presented. 

We expect that in an NNLO (^-expansion calculation the Delta's effect on f3 p will roughly 
double. a p will also be modified, thanks to the graphs in Fig. 4 and the w-channel Delta- 
pole graph with two qe couplings. Estimating these effects gives the theory error bars which 
appear in the second line of the table. We note that, even though NNLO effects make a 
significant difference in the values for a p and fl p , their impact on the low-energy differential 
cross section is not large, being represented by the theoretical error band in Fig. 7. This 
suggests that the extraction of polarizabilities from Compton data is a delicate process. 



Reference 




Pp 


NLO HB^PT [19] 
NLO 5 [this work] 
NLO SSE [25] ([23]) 


12.2 
10.21" 
16.4 (20.8) 


1.2 

9.1 (14.7) 


PDG average [39] 
LEGS [4] 
MAMI [5] 
Beane et al. [20] 


12.0 ±0.7 

11.8 ±2.0 

11.9 ±2.1 

12.1 ± 1.6 


1.6 ±0.6 
1.4 ± 1.5 
1.2 ± 1.4 
3.2 ± 1.2 



TABLE II: Proton polarizabilities in xPT, the <5-expansion, and the SSE, compared to values 
extracted from experiment. Results are in units of 10~ 4 fm 3 . 



VI. CONCLUDING REMARKS 

The expansion developed in this paper for the 7N amplitude is based on the scale hierarchy 
to^ < A < A. The EFT expansion parameter employed here is 5, which represents both 
the ratio of m n to A and the ratio of A to A. The success of the resulting "^-expansion" is 
to be judged by its efficacy as an EFT of Compton scattering. 

The existence of two different low-energy scales in the theory forces us to develop in- 
dependent power countings for the two different photon-energy regimes. For uj ~ m n the 
Compton amplitude obtained is exactly that of up to effects from the Delta which are 
down by 5 3 relative to leading. On the other hand, if uj ~ A, OAR graphs dominate the 
Compton amplitude. Resummation of the Delta self-energy is necessary in these graphs. 
The Delta propagator appearing in them acquires a finite, energy-dependent, width. Up to 
effects suppressed by S 3 this is the only correction to the Delta propagator that is necessary. 
Vertex corrections also appear at NLO in this region. 
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We performed a calculation of the 7N amplitude which includes all effects that are of 
leading or next-to-leading order in the region uo ~ m n , as well as all effects that are of 
leading or next-to- leading order in the region u ~ A. The sum of all of these pieces of 
the amplitude defines our NLO calculation. Thus, in each of the two regions considered, 
mechanisms of NNLO in that particular region are included in the calculation. Nevertheless, 
effects of relative order <5 3 (e.g. A-7T loops) are omitted in the low-energy region, while effects 
of relative order 5 2 (e.g. two-loop dressing of the Delta) are omitted in the higher-energy 
domain. These neglected effects define the theoretical error bar of our calculation. 

Note that in this calculation tree-level graphs are computed relativistically, while for loop 
graphs we use the HB^PT result. While this gives the correct 7N amplitude up to the order 
to which we work, a fully-relativistic treatment of 7rN loops would be more aesthetically 
pleasing. 

After fitting the only two free parameters in our theory, the E2 and Ml 7./VA couplings, 
good agreement with the low-energy 7p data is found. The spin-independent polarizabilities 
a and f3 that result from our NLO calculation are in reasonable agreement with contemporary 
extractions from data [5, 20]. 

The development of the (^-expansion opens up a number of avenues for further study. 
Higher-order calculations of Compton scattering on the nucleon will be necessary to see if the 
good agreement found at NLO persists, and if the expansion in powers of S is well-behaved 
or not. Also, the use of the (^-expansion in other processes, e.g. pion photoproduction, is an 
important potential future application. Indeed, the EFT presented here should ultimately 
be judged by its success in simultaneously describing data on nucleon Compton scattering, 
pion photoproduction, and 7rN scattering. Only then can the reliability of the EFT, and 
hence of our extraction of g& and qm from 7N data to NLO in the (^-expansion, really be 
judged. 

The use of the (^-expansion in two-nucleon systems might also reap significant rewards. 
We plan to use the amplitude developed here, together with consistent two-body currents, 
in a calculation of 7D scattering [42]. It is also possible that the 5 expansion could be prof- 
itably employed to organize the Delta contributions to the chiral nucleon-nucleon potential 
developed in Refs. [43-45]. 
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APPENDIX A: RESULTS FOR COMPTON INVARIANT AMPLITUDES 

We define the photon energy as: u = (p • q)/M N = (s - M%)/2M N . The results for the 
invariant amplitudes are presented in the nucleon mass units (M N = 1). Below Z — 1 for 
the proton, Z = for the neutron. Expressions for the Delta-exchange are obtained using 
the algebraic manipulation program Form [46]. 
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1. Nucleon s-channel: 



Ai(u,t) 

A 2 (u,t) 
A 3 (u>,t) 
A 4 (u;,t) 

A 5 (w,t) 
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A s (u,t) 
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3. Delta s-channel: 
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4>(«;,i) = F(uj,A){1(2 + A) G 2 M -luj [2(1 + A) G 2 M + AG 2 E + (2 + A) G M G E 



21 



— lo 2 Gm (Gm — Ge) + \ t GmGe + \ A Gm (Gm + Ge) — w Gm {Gm — Ge) } 
A 6 (lu, t) = F(u>, A) {lu G m [G m + \ A (G M + G E )\ + cu 2 G M (G M - G E ) 

+ \t G 2 M + 5 A Gm (Gm + G#) + to Gm (Gm — G\e)] j 
A 7 (w,i) = F(o;,A){(2 + A)(^Gm + ^G m G b )-^AG| 

+ to G 2 M ~~ 12 G\ ~~ \ GmGe} + jq t (Gm + G^) 2 j 
A 8 (u,t) = F(u, A) {|(2 + A) G 2 M + cu [4(1 + I A)G M - | (2 + A) G M G E 

+ \u 2 (G 2 M + G 2 E-AG M GE)-\tG 2 M(\ + uj)} 



with G m ,e = 2 (2 3 +A) 9m,e and 
F(c,A) = - 2 



3 s - Ml + iM A r(s) 32u- A(2 + A) +iM A r(u) 



(A4) 



where f = T 3 t T 3 is the isospin factor, and the width r is given by Eq. (42). 

The corresponding w-channel graphs are obtained by crossing so that the crossing- 
symmetric amplitude is given by 



A { (uj,t) + Ai(-u',t), fori = 1,2, 8 
Ai(u,t) — Ai(—u' ,t), for i = 3, . . . , 7 



with uj' — p ■ q' — (1 — u)/2. 



4. HBxPT vriV loops 



Defining C = (^fs-M N )/m n and r = -2( 2 (1 — cos#), where 9 is the center-of-mass angle 
between the incoming and outgoing photon momenta, one finds [19]: 



A 1 (s,t) 
A 2 (s,t) 
A 3 (s,t) 
A 4 (s,t) 
A 5 (s,t) 
A 6 (s,t) 



9 2 A m « L f, Ji, 



2-r 



'— r 



1 v^r 
-arctan — A(C^) 



2-r 



87r/ 2 m^ (-r)3/2 



[/i(C,r)-/ 2 (C,r)], 



8tt 2 / 2 
9l 



47T 2 /2 m7r 

8tt 2 />. 



-arcsin 2 C-C + 2C 4 sin 2 £/ 3 (C,r) 
h(C,r), 



[/ 5 (C,r)-2C 2 cos^ 3 (C,r)], 



87r 2 /2 m?r 



[J 5 (C,r)-2C 2 /3(C,r)], 



(A5) 
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where 



with 



h(C,r) 

HC,r) 
h(C,r) 
h((,t) 



1 , (l-^v 7 ^ 1 
az arctan — , = 

2y/l - ( 2 Z 2 

1 2(1-Z)y/-T(l-<?Z*) 



= L dz Mi 



4(1 - ( 2 z 2 ) - t(1 - z) 2 ' 

i f 1 1 x(l -x)z(l -zf 

ax / az — 

o jo o d 

1 (fa / 1 (i ,£<i_i) arcsi 4, 

o Jo S R 

1 , f l .{l-zf . (z 

ax / az arcsm—, 

o jo S R 



(z (zS' 
arCSm i? + ^ 



S= Jl- ( 2 z 2 - r(l - z) 2 x(l -x), R=Jl~ r(l - ^) 2 a;(l - x). 



(A6) 
(A7) 



APPENDIX B: RESULTS FOR POLARIZABILITIES 



The nucleoli electric {oln) and magnetic (/3at) polarizabilities: 



a N = 



®j^ A (NB) (00) A (NB) (00) _^(jA_\ 2ag 2 E 



^ = -4^(0, 0) = 2^r (0,0) ' (M , + M A)2 A 



7ra / p A 



+ 



2«Pm 



, (Bl) 



where Aj 7VB ' ) are the amplitudes with the Born graphs subtracted. 
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FIG. 7: Angular distribution of the 7p differential cross section at low energies. Data points are 
from Ref. [5]. The long-dashed orange line represents the sum of nucleon and pion Born graphs, 
the blue dashed line gives the NLO xPT prediction, and the red solid line is the full result at NLO 
in the (5-expansion. The dots give an estimate of the theoretical error. 
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FIG. 8: Angular distribution of the 7p differential cross section. Legend for the curves is the same 
as in Fig. 7. Data points are from Refs. [5] - MAMI'01, [1] - SAL'93, [4] - LEGS'01. 
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FIG. 9: Energy dependence of the 7p differential cross section at 90°. Data points are from Refs. [1] 
- SAL'93, [3] - LEGS'97, and [4] - LEGS'01. The red solid, orange long-dashed, and blue dashed 
lines are as above, while the purple curve now gives the contribution arising solely from the Delta, 
and the green short-dashed curve is the effect of 7rN loops alone. 
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